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Abstract 

The dynamical heterogeneities (DH) in non-ergodic states of an attractive colloidal 
glass are studied, as a function of the waiting time. Whereas the fluid states close 
to vitrification showed strong DH, the distribution of squared displacements of the 
glassy states studied here only present a tail of particles with increased mobility 
for the lower attraction strength at short waiting times. These particles are in the 
surface of the percolating cluster that comprises all of the particles, reminiscent 
of the fastest particles in the fluid. The quench deeper into the attractive glass is 
dynamically more homogeneous, in agreement with repulsive glasses (i.e. Lennard- 
Jones glass). 
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1 Introduction 

Two distinct glasses have been predicted and identified in hard spheres with 
short range attractions: a repulsion driven glass at high density, and an at- 
traction driven one, formed at low temperatures (or high attraction strength) 
[1,2]. Whereas the former one is caused by the steric hindrance of the particle 
cores and the so-called cage effect, the latter forms due to the bonding be- 
tween particles. This system is realized experimentally by a colloid-polymer 
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mixture, where the effect of the polymers is to induce an effective attraction 
between the colloids [3]. Both glasses have been indeed identified, although 
the attractive one, which at low concentrations is termed 'gel', often competes 
with (and inhibits) liquid-gas phase separation [4]. 

Dynamical heterogeneities (DH) have been found in the proximity of repulsion 
driven glass transitions by computer simulations, i.e. in Lennard- Jones mix- 
tures [5,6], or hard spheres [8]. In these cases, while the system is structurally 
liquid-like (homogeneous), a population of particles of increased mobility is 
observed. As the glass transition is approached from the fluid side, the het- 
erogeneities become more pronounced, but decrease again deeper in the glass 
[7,9,10]. The role of these dynamical heterogeneities in the glass transition is 
as yet unclear; whereas mode coupling theory focusses on averaged quantities 
and neglects them [11], the so-called facilitated dynamics theories give DH the 
central role for their description of the glass transition [12]. 

In recent works, it has been shown that DH can be found also in attractive 
glasses, by studying the distribution of particle displacements in the system 
[13,14,15]. In fluid states close to the transition two populations of parti- 
cles were found, separated by a minimum in the displacement distribution. 
A similar feature has been found also in repulsive glasses, which could imply 
a common origin [6,15]. However, the low density of the attractive glass, as 
low as 4> c = 0.4, causes structural heterogeneities as well; the system forms a 
percolating cluster of high density material, leaving voids with no particles. A 
correlation between structural and dynamical heterogeneities is thus possible, 
showing that 'fast' particles are in the surface of the cluster, whereas the 'slow' 
ones are mostly trapped in the inner parts of it [13]. 

In this work, we study the DH inside the non-ergodic region, for two differ- 
ent states, and compare them with those of the equilibrium systems. Only 
one population of particles can be identified from the distribution of particle 
displacements, and the distribution is narrower for the state with stronger 
attractions. Moreover, as the systems age, they become more and more homo- 
geneous, from the point of view of the dynamics. Both results indicate that 
the strongest DH are obtained in the fluid side of the phase diagram, close to 
the glass transition. As a side remark, it must be noted that the structural 
heterogeneities mentioned above persist in the out-of-equilibrium systems, and 
thus are not the sole origin of the DH in attractive glasses. 



2 Simulation details 

We have performed computer simulations of a system composed of 1000 soft 
core (V c ~ r~ 36 ) particles with attractive interactions given by the Asakura- 
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Oosawa (AO) potential [3]. It models a mixture of colloids with non-adsorbing 
polymers, and the range of attraction is set by the polymer size. In order to 
guarantee full access to the whole parameter space, phase separations have 
been inhibited. Crystallization is avoided by polydispersity (flat distribution, 
5 = 10% width), and liquid-gas demixing by a repulsive barrier extending to 
two mean diameters. Further details of the interaction potential can be found 
in previous works [14]. 

Length is measured in units of the average radius, a, and time in units of 
^4a 2 /3v 2 , where the thermal velocity, v, was set to ^4/3. Equations of motion 
were integrated using the velo city- Ver let algorithm, in the canonical ensemble 
(constant NTV), to mimic the colloidal dynamics, with a time step equal 
to 0.0025. Every n t = 100 time steps, the velocity of the particles was re- 
scaled to assure constant temperature. The range of the attraction is set to 
0.2a. Density is reported as volume fraction, C = 4/3ira 3 (l + (5/a) 2 ) n c , 
with n c the number density, and the attraction strength is measured in units 
of the polymer volume fraction <p p (at contact the AO attraction strength is 
16kBT<pp). The attractive glass transition for this system has been studied 
previously [14,16]. An MCT analysis of the results (diffusion coefficient, time 
scale and viscosity) yields a transition point at 4> p = 0.4265 for the colloid 
density <p c = 0.40. 

For the study of aging here, the systems were equilibrated without attraction 
((j) p = 0) at C = 0.40, and then instantaneously quenched to the desired (f> p 
at zero time, t — 0. Two attraction strengths have been studied, <p p = 0.50 
and P = 0.80, lying beyond the nonergodicity transition. In both cases, 25 
independent simulations have been performed, and the evolution of the system 
has been followed as a function of the time elapsed since the quench, called 
waiting time, t w . Correlation functions thus depend on two times: t w and 



3 Results and discussion 

In fluid states close to the attractive glass, increasing DH have been found, 
the stronger the attraction [13]. Two populations of particles are observed to 
appear as the attraction is increased, one of mobile particles and another one 
of quasi-immobile particles (see Fig. 1) [14]. (The minimum in the distribution 
of squared displacements allows for an unambiguous identification of almost 
every particle.) These two populations are structurally well differentiated: the 
particles in the 'skin' of the gel, with a small number of bonds, are mobile, and 
the particles in the inner parts of the gel, the 'skeleton', are quasi-immobile. 
The populations were observed to be stable for long times, although particles 
can change from one to the other. It is thus an equilibrium feature of the 
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Fig. 1. Distribution of squared displacements at different times: from left to right, 
t = 100, 200, 400, 800, and for different polymer fractions, cp p , as labeled. Note 
that as 4>p increases, two populations of particles with different mobilities appear 
in the system. The glass transition is located at <p p = 0.4265, estimated from MCT 
analysis (power law fittings) [14]. 
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Fig. 2. Mean squared displacements for 4> p = 0.50 and 4> p = 0.80, as labeled, for 
different waiting times; from left to right, t w = 8, 16, 32, 64, . . ., 16384. The dashed 
line shows the MSD for hard spheres at this density. 

system, partly induced by the structural heterogeneity. 



All of the states presented in Fig. 1 are, nevertheless, fluid states. The struc- 
tural properties do not depend on time, nor the dynamical ones on the initial 
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Fig. 3. Evolution of the mean number of neighbours per particle in the system after 
the quench for both states, as labeled. 




Fig. 4. Snapshot of the system with 4> p = 0.50 at a long time, t = 30000. (Qualita- 
tively similar snapshots are obtained for <j) p = 0.80.) 

time, and the correlation functions averaged throughout the whole system de- 
cay to zero. By further increasing the attraction strength, however, the system 
falls out of equilibrium and shows aging. In Fig. 2, the mean squared displace- 
ment (MSD) is presented for the two quenches studied here, (f> p = 0.50 and 
4> p = 0.80; different lines are the MSD for different waiting times, t w . As t w in- 
creases, bonds are formed between particles, as shown in Fig. 3 for both states, 
which hinders the motion of the particles, causing the dynamical arrest. 
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Fig. 5. Distribution of squared displacements for (p p = 0.50 and 4> p = 0.80 and three 
waiting times, as labeled. The different lines represent the distribution at t = 100, 
200, 400, 800, 1600, 3200, 6400, 12800 and 25600, from left to right, respectively. 



In the MSD, a plateau develops at short distances, signaling the localization 
length, clearer in the case of the higher p , where the bonds are stronger and 
the cage (network of bonds) of the particles is tighter, as shown by the shorter 
localization length (implying higher non-ergodicity parameters in the density 
correlation function). The mean number of bonds per particle is, however, 
lower. Also, the aging is more dramatic and at long waiting times the MSD 
hardly reaches the range of the attraction, contrary to previous works where 
the plateau in MSD, or the density correlation function, was not observed [17]. 
Here we will concentrate on the DH, and a full analysis of the aging will be 
presented elsewhere [18]. 

A snapshot of the system with P = 0.50 is presented in Fig. 4. The system 
forms an intricate structure, with voids and tunnels, similar to the fluid states 
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Fig. 6. Squared displacements the particles during time t, as a function of its mean 
number of neighbours during this time, for (j) p = 0.50 and t w = 1024. (Only ten 
simulations are considered, i.e. 10000 particles). Note the vertical bands at integer 
numbers of neighbours (especially for t = 1600), due to particles that have not 
changed their neighbours. 



presented above. Therefore, the "skeleton and skin" picture presented above, 
with two different populations of particles, could still be applicable. In Fig. 5 
we present the distribution of squared displacements at different times for both 
quenches. Strikingly, the distribution is monomodal in all cases, and the peak 
evolves to larger displacements as time proceeds. The strong DH observed 
in the fluid states, have, therefore, disappeared in the glass side. Moreover, 
the deeper quench, <p p = 0.80, shows narrower distributions at all waiting 
times than the quench at <f) p = 0.50, indicating that the dynamics is more 
homogeneous the deeper the state is in the glass, in agreement with findings 
for Lennard- Jones (repulsive) glasses [10]. 

At <f) p = 0.50, however, a tail in the distribution to long distances can be ob- 
served, caused by some particles that can travel long distances; the number 
of which decreases with t w . This feature is reminiscent of the population of 
fast particles observed in the fluid (Fig. 1). The origin of this population is 
studied in Fig. 6, where the squared displacement of every particle for a given 
time is correlated with the mean number of neighbours of the particle during 
this time. The plot shows that indeed the fastest particles in the system have 
less neighbours on average, whereas the particles with many neighbours, move 
very little. Therefore, the simple picture of "skeleton and skin", can still be 
applied for the attractive glass close to the transition. The number of these 
fast particles is, nevertheless, decreasing with waiting time, as observed by 
comparing similar times t for different waiting times t w in Fig. 5. Accordingly, 
the mean number of bonds per particle increases (see Fig. 3), implying a com- 
paction in the system, and thus having fewer particles in the skin, and even 
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those, more tightly trapped. However, we cannot state whether the population 
of fast particles will reach a steady state or if it will vanish eventually at very 
long waiting times. At <p p = 0.80, on the other hand, the tail of fast parti- 
cles is absent in the distribution of squared displacements, although similar 
structural heterogeneities are observed. 



4 Conclusions 

We have shown that the strong DH found in fluid states close to the attrac- 
tive glass transition in colloids with short range attractions decrease again 
deep into the non-ergodic region. The distribution of squared displacements is 
monomodal and no particles with increased mobility are observed (0 P = 0.80). 
However, close to the glass transition, in the glass side, aging is slower and 
some DH can still be detected: a tail in the distribution of squared displace- 
ments indicates fast particles, that can be identified with particles in the outer 
parts of the particle network. This feature is reminiscent of the population of 
fast particles found in the fluid states close to the transition. However, the sta- 
bility of this population of fast particles in the glass cannot be established, but 
should not be present in a truly arrested glassy state. The results presented 
here agree with those of repulsive glasses, contrary to other comparions be- 
tween attractive and repulsive glasses [17]. 
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